Logarithmic periodicities in the bifurcations of type-I intermittent chaos 
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The critical relations for statistical properties on saddle-node bifurcations are shown to display 
undulating fine structure, in addition to their known smooth dependence on the control parameter. 
A piecewise linear map with the type-I intermittency is studied and a log-periodic dependence is 
numerically obtained for the average time between laminar events, the Lyapunov exponent and 
attractor moments. The origin of the oscillations is built in the natural probabilistic measure of the 
map and can be traced back to the existence of logarithmically distributed discrete values of the 
control parameter giving Markov partition. Reinjection and noise effect dependences are discussed 
and indications are given on how the oscillations are potentially applicable to complement predictions 
made with the usual critical exponents, taken from data in critical phenomena. 

PACS numbers: 05.45.Pq, 05.45.-a, 64.60.Fr 



Dynamical bifurcations, when the qualitative behav- 
ior of a natural phenomena changes due to the variation 
of a control parameter, is a subject described by a well 
founded mathematical theory y|. Their extension from 
equilibrium thermodinamic phase transitions to the bi- 
furcations in non-equilibrium systems is applied within 
and outside physics 0,0]. Some bifurcations are abrupt, 
like first order equilibrium transitions, and others fol- 
low critical relations with a continuous dependence in 
the control parameter. Sufficiently close to the bifur- 
cation the system properties follow characteristic func- 
tions giving rise to critical exponents which are signa- 
tures of the phenomena. Log-periodic oscillations modu- 
lating critical functional relations have been reported in 
studies of bifurcations 0- Q| ■ As reviewed by Sornette 

0, many ongoing research on applications are active, 
among them the prediction of catastrophic events. Tran- 
sient chaos p| and other escape phenomena |(| also dis- 
play such log-periodicities. In simple maps, the period 
doubling Feigenbaum cascade is an early example where 
log-periodic behavior appears .7]. In non attracting sets 
of two dimensional maps, the topological entropy have 
been shown to present log-periodic oscillations 8]. 

Among the dynamical bifurcation phenomena many 
can be cast on the class of intermittent chaos, ranging 
from the onset of turbulence to synchronism of chaotic 
systems [Tol J T 3 . As originally proposed by Pomeau and 
Manneville @>"E1 113 > these instabilities can be modeled 
with simple one-dimensional maps and classified as type- 

1, II, and III, according to the Floquet multipliers of the 
map crossing the unity circle in the complex plane at 
1, at a pair complex conjugate values or at —1, respec- 
tively. Type I, which will be discussed here, occurs when 
a saddle-node or tangent bifurcation is approached. It- 
erations of the map near the value of the virtual fixed 
points are identified with the laminar events while the 
reinjection iterates correspond to the turbulent bursts 
that occur in a erratic manner. Such bifurcation has 
been reported in many experimental systems |l'il Il4j 
and it is abundant in the logistic and other simple maps. 
Renormalization group procedure to deal with type-I in- 



termittency has already been used 0] and give analyti- 
cal asymptotic expressions for statistical averages in the 
maps. These are treatments that only give smooth criti- 
cal dependences. 



The tangent bifurcation has analogy with second or- 
der thermodynamic phase transitions. The continuous 
scale invariance of thermodynamical phase transitions 
gives smooth monotonic critical relations for the statis- 
tical thermodynamical variables or the susceptibilities, 
as the order parameter approaches the transition. How- 
ever, as mentioned above, oscillations do occur in the 
corresponding properties of non-equilibrium dynamical 
bifurcations. Discrete scale invariance in the dynamical 
systems bifurcations underly the oscillating features 0. 
In type-I intermittent chaos with the tangency having 
second and fourth order nonlinearity, they appear with 
power-law periodicity 0, , instead of log-periodicity. 



Herein, to study the origin and properties of the log- 
periodic undulations, the simplest unidimensional map 
exhibiting chaos and a saddle-node bifurcation was con- 
sidered. It is a piecewise linear map first proposed to 
study critical scaling laws of type-I intermittency by 
Shobu et al. ^^|. A reason to take this map is the exis- 
tence of an analytical approach to find the discrete values 
of the control parameter at which the scaled shape of the 
map function reproduces itself. The average length of 
laminar phases, the spectrum of fractal dimensions, the 
Lyapunov exponent, and the statistical moments of the 
chaotic variable, when calculated numerically, show log- 
periodic oscillation. Such numerical results correspond to 
time averages and became possible due to the great com- 
putational power achieved with current digital machines. 
They can be compared to, and surpass in detail, analyt- 
ical results from averages in ensembles valid at specific 
values of the control parameter. The oscillation reported 
here are one such case. Unfortunately we were unable to 
extend the theory of renormalization group |15| to pre- 
dict such periodic property. 
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The equation for the map is E3 



ax n + b, 



for x„ < 



x n+l — 



l [x n - (1 - 6)] + 1, for C < X n < 1 - b 



-b-\x n - 1), 



for 



> 1-6 



(1) 

where < b < 1/2 is a constant , a = 1 — 26 + e > 0, 
and c = (1 — 6)/(l + a) > are expressed as functions 
of the control parameter e. The map is represented in 
figure [2 with b = 0.2 and e = 0.0332. When e < the 



n+1 




FIG. 1: The So-Ose-Mori map with an unstable periodic or- 
bit, calculated with b = 0.2 and e = 0.0332. 



first two linear segments cross the diagonal and a sta- 
ble fixed point exists at the first crossing. At e = the 
corner point of the first two linear segments touch the 
diagonal, at x = 1/2, and a bifurcation occurs similar 
to a saddlc-nodc when two fixed points disappear. For 
e > 0, only the third linear segment intersects the diag- 
onal with an absolute value of slope greater than 1 and 
the map is chaotic. Iterates at left of (1 — 6) shall be 
considered as laminar, and iterates on the negative slope 
will be referred to as reinjection events. With e > 0, c 
is approximately proportional to the minimum distance 
between the map and the y — x diagonal. 
Shobu et al. 

Hi 

have shown that, for this map, sim- 
ple Markov partitions exist at discrete values of e. When 
this happens, as one iterates the map, the probability of 
visits, or the natural measure, over the unit interval is a 
countable set of step intervals with constant value. The 
map dynamics is that of a finite Markov chain, diverging 
for e — > 0, and it is analytically solvable. The appearance 
of unstable periodic orbits, reinjecting at x = and it- 
erating 2n + 2 times until reaching precisely the value 1, 
mark the condition for the Markov partition with 2n + 1 
constant steps for the natural measure. Figure Q shows 
the map and the iterates of such an orbit with period 12. 
The specific values of e = e n for which the orbit of period 
2n+2 exists follow the recursive relation e„ = (1— 26)e ra _i 
|l8j . and, as good approximation for n > 3, 



'n=-2b[i- 7 ^)(l-2b) n - 1 



(2) 



The factor (1 — 26) will appear as the period of the fine 
structure in logarithmic scale. At the values of e n one can 
obtain analytically the ensemble averages for the average 
length of laminar events, (Z), for the first order moment of 
the attractor, (x), |l8j. and for the Lyapunov exponent, 
(A) , The analytical expressions first determined by Shobu 
et al. JL8], and resulting from the summation with the 
stepwise uniform probabilistic measure, are 



<Z) = 



1 hfl 
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The ensemble averages must coincide with time aver- 
ages Z, A and x, calculated using iterations of the map and 
starting from typical initial conditions. The length of a 
laminar event is considered as the number of iterations 
"in the channel" before a reinjection. Any consideration 
of a narrower range around the corner point as the chan- 
nel causes just a shift in the value of I for a given e, but 
the slope and fine structure dependence on e remains the 
same. The Lyapunov exponent (we use A for A and(A}) 
was calculated numerically using a large number of iter- 
ations, N, with the expression 



1 N 

A = ^£H./>n)l 



and the average of the dynamical variable x with 



1 r 

x = n ^ X% 

71=1 



(6) 



(7) 



Figure shows the average length of laminar events 
calculated both with time averages and ensemble aver- 
age. In the last case, only a discrete set of values of e 
from eq.|21can be done. The values of Z decrease with the 
control parameter e. Thus if one calculates and inter- 
polates the analytical discrete values E3 one never gets 
the periodic effects. Conversely, for the time averages the 
values of e can be picked at will. In the figure for Z, 2000 
values of the control parameter e were equally spaced on 
a logarithmic scale and at each one, N — 3 x 10 7 itera- 
tions were done for the averages, with double precision 
floating point arithmetics. One observes clearly that, in 
addition to the smooth power-law divergence as hie -1 , Z 
is decorated with a modulation whose period is given by 
the positions at which each Markov partition is formed. 

Numerical computations of the average length of lam- 
inar events and average Lyapunov exponent are show in 
fig. |3J Logarithmic scale is used in the absciss and the 
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FIG. 2: Average length of laminar events. The dots are an- 
alytical ensemble averages and the continuous curves are nu- 
merical time averages, iterating the map, as discussed in text. 




FIG. 3: Average length of laminar events (a), and (b) aver- 
age Lyapunov exponent, A, in the SOM map. The dots are 
analytical ensemble averages and the continuous curves are 
numerical time averages. 



critical behavior of I has a linear smooth envelope behav- 
ior as hie -1 for small e. Furthermore, its oscillations are 
periodic as function of lne. The behavior for (x) is pro- 
portional to the one exhibited by A, as one expects com- 
paring eqs.0]and[S] Changing b on equation ^ changes 
the spacing between successive e n , and the correspond- 
ing period obtained in time averages follows consistently 
with the scaling factor e„/e„ + i = (1 — 26). Therefore, nu- 
merical artifacts for log-periodicity are completely ruled 
out in the computations. This b dependency is shown in 
fig. 01 The same oscillating features were obtained (not 
show here) for D q (e), the fractal dimension of the attrac- 
tor, calculated for q = 0, 1 and 2. The critical exponent 
of the smooth envelope of any of the averages do not 
change with b. 

Calculations were also done using a random reinjection 
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FIG. 4: Average of the dynamical variable in the SOM map 
with different scaling factor: (a) b = 0.2 and (b) b = 0.3. 



value in a map with the two first linear segments as the 
laminar phase part of eq. ^ As long as the reinjection 
is non singular, oscillations with the same period are al- 
ways obtained. This fact is indication of the local nature 
of the oscillations. While the Markov partitions of the 
SOM maps are a global result, their relevant contribu- 
tions for statistical averages near criticality come from 
the uniform density steps near the narrow channel. Any 
approximately uniform reinjection on this region suffices 
to give oscillating averages. 

Finally some discussion on noise. The effects of noise in 
type-I intermittency have been studied analytically [l5| 
and experimentally |l9j |. Here, a term of additive gaus- 
sian white noise with amplitude D was added to the right 
hand side of eq. ^ Fig. shows 7, when the noise ampli- 
tude was D = 10 -6 . For e 3> D one expects the dynam- 




FIG. 5: Averages length of laminar events in the SOM map 
with additive noise. 

ics to remain unaffected. When the separation between 
successive e„ becomes of the order of the noise ampli- 
tude, D, the corresponding and smaller oscillations are 
blurred out. For even smaller e (< D), the number of 
iterations in the narrowest part of the channel no longer 
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depends on e, because a single noise event may cause a 
jump across the waist of the channel, which, otherwise 
would have a number of iteration proportional to lne -1 . 
Therefore, I saturates. Such noise sensitivity is shown 
for x in fig. |B| It is relevant to observe that, consistent 
with the above explanation, in both figures, and the 
oscillations cease to exist for e one order of magnitude 
higher than the value at which the slope of the average 
saturates. One also may see a smoothing of the ultra- 
fine structures on the undulations. Smoothing, erasure 
of the undulations, and saturation of the amplitude en- 
velope also occurs in the Lyapunov exponent when noise 
is present. It may also be observed in the figures that 
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FIG. 6: Average of the dynamical variable, x, in the SOM 
map with additive noise. 

the amplitude of the undulations on I and x are bigger 
at different ranges of e. This might be used to justify 
the importance of acquisition of different data series and 
different averaged quantities in experiments. 

To conclude it must be stressed that oscillations in 



the critical behavior of statistical properties occur univer- 
sally in the tan gent bifurcations of unidimcnsional nor- 
mal form maps ]l3| and in the logistic map 16]. They 
also occur in the saddle-node bifurcations of continuous 
chaotic systems, as we verified numerically integrating 
the Rossler, the Lorenz and the equations for a laser with 
saturable absorber [2(j. However none of these follow 
the log-periodic dependence. The map herein was cho- 
sen to stress the relationship between log-periodic fine 
structures, existence of Markov partition and self similar 
scalling. From a mathematical point of view much more 
remains to be done on this connection. The undulations 
are important features of bifurcations that can be subject 
to experimental tests and applications. To the present, 
in experiments with intermittent chaos only the smooth 
critical behavior have been studied. The observation of 
log-periodic oscillations has been proposed as a technique 
to anticipate the criticality in natural phenomena, partic- 
ularly in earthquakes 3] . The simple map presented here 
is far from modeling earthquakes but it is an workable 
example to show how the critical value of a parameter 
controlling bifurcations might be predicted ahead of crit- 
icality, by complementary fitting of oscillations in the fine 
structure at higher values of the control parameter, along 
with the usual fitting of the monotonic critical behavior. 
It also shows how improvement in extracting informa- 
tion on a bifurcation with intermittent chaotic should be 
enhanced by measuring different statistical variables on 
different ranges of the control parameter. 
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